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1. INTRODUCTION 

Multilevel inverters (MLIs) are a popular power electronics device used in mid-voltage and 
high-power applications instead of conventional two-level inverters. Basically, there are three conventional 
topologies for MLIs these are; cascaded H-bridge MLI (CHB-MLI), capacitor-clamped MLI (CC-MLI), and 
diode clamped MLI (DC-MLD [1]-[3]. The robust and modular structure of CHB-MLI made it the most 
commonly used topology among the existing topologies [4]. Conventional inverters have high efficiency and 
a high-power ratio. However, the semiconductor switches experience higher stress as the desired power and 
voltage value are increased. Today, it is not possible to produce semiconductor power switches that can 
withstand medium voltage levels. Therefore, the use of MLI with lower voltage stress in switches in 
medium/high voltage applications has become popular [5]-[8]. 

However, this technology faces a harmonic problem that causes harmful effects on the load. Among 
harmonics, the most dangerous harmonics for the system are low-order harmonics because these harmonics 
cause torque fluctuations in electric motors, decrease efficiency, and reduce the lifetime of the system [9]-[12]. 
The control and modulation methods used by these inverters have become more challenging as MLI topologies 
have advanced [12]-[14]. An essential issue that needs to be resolved is the existence of low-order dominant 
harmonics in the cascading output voltage generated by MLI. These harmonics mostly cause malfunction, 
increased loss, and voltage ripple. These harmonics also have an impact on power quality [15]. Despite the fact 
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that there are other controls and modulation techniques, the selective harmonic elimination (SHE) technique is 
the primary modulation technique that solves the aforementioned issues. 

Hence, the solution to the problem is extremely difficult. Three sorts of techniques can be utilized to 
address SHE issues; these are numerical techniques, algebraic techniques, and meta-heuristic algorithms. The 
gradient optimization and Newton-Raphson (NR) [16] methods can be given as examples of numerical 
methods; however, these iterative methods are mainly sensitive to the initial value and divergence problems are 
likely to occur, especially as the inverter level increases. If the initial value is not selected properly, the number 
of iterations may increase greatly or the solution may not be found. Algebraic methods are techniques such as 
resultant theory and Groebner bases theory for calculating optimized switching angles. These techniques 
involve the conversion of transcendental equations into the corresponding polynomial equations. These methods 
do not depend on the initial value but are very complex computationally, therefore, these methods can be 
applied only to low-level inverters. In meta-heuristic methods, SHE equations are solved using nature-inspired 
methods such as genetic algorithm (GA) [17] particle swarm optimization (PSO) [17], [18], LSHADE/EpSin 
technique (LSHADE) [19], whale optimization technique (WOA) [20], [21], and enhanced krill herd (EKH) 
[22]. There are certain common problems of algorithms such as PSO and GA; these include convergence to 
local optimum, early convergence, and either only global solution or only local solution. 

This study proposes a novel red deer algorithm (RDA)-based optimization approach to the SHE 
problem that can overcome the above problems. The proposed RDA optimization has been applied for the 
first time for the SHE problem [23], [24]. According to the optimization methods used for the SHE problem 
in the literature, the most obvious superiority of RDA optimization is that it can adjust the diversification and 
intensification phases. Owing to this feature, RDA can trade off between the optimal minimization total 
harmonic distortion (OMTHD) and SHE (it can use the local and global search together). To validate the 
concept, switching angles for OMTHD and SHE-pulsewidth modulation (SHE-PWM) using RDA are found 
and applied to a 3-phase 11-level CHB-MLI. In addition, the harmonic results of the RDA-based SHE 
technique have been compared with the results of other techniques in the literature. 


2. CASCADED H-BRIDGE MULTILEVEL INVERTER 

Figure 1(a) displays the circuit diagram for the N-level three-phase CHB-MLI. Each phase is made 
up of a single-phase H-bridge inverter that is connected in series and is powered by a separate, discrete DC 
source. If the voltage value of DC sources is equal, it is called a symmetrical MLI, otherwise, it is called an 
asymmetric MLI. The proposed algorithm in this study is applied to the symmetrical structure. Figure 1(b) 
depicts the AC output voltage waveform (Vo) for an N-level MLI supplied by a symmetrical source. The 
output voltage for each phase can be expressed as (1): 


V = Voc (1) 


Where k is the number of discrete DC sources. Also, k represents the number of H-bridge modules. The 
relationship between the number of DC sources (k) and the voltage level N can be given as (2): 


N=2k+1 (2) 
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Figure 1. Three-phase N-level CHB-MLI (a) circuit diagram and (b) output voltage waveform for single 
phase 
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3. SELECTIVE HARMONIC ELIMINATION PROBLEM IN MULTILEVEL INVERTERS 

Low-order harmonics are eliminated by the SHE-PWM switching method via the creation of 
switching angles at the fundamental frequency. The number of discrete DC sources in this approach equals 
the number of switching angles (k). Switching each module just once will eliminate k-1 maximum number of 
harmonics from the output voltage of the inverter. The output voltage can be fourier expanded to reflect the 
required nonlinear harmonic equations to attain the optimal switching angles. The expression of the output 
voltage that contains all the harmonic components is given as (3) [5], [6]: 


Vo (wt) = Elok n=1,3 toc Tos (n01) + cos (n83) + = + cos (n8y)]sin (wt) (3) 


ite nn 


For the output voltage, the harmonic components are mostly composed of the fundamental 
frequency and its multiple. The symmetrical output voltage represented in Figure 1(b) ensures the 
non-availability of even harmonic components. The output voltage consists of only odd harmonic 
components as expressed in (3), where w=2rf and f is the fundamental frequency. The expression of the 
fundamental harmonic component (V1) of the output voltage is given thus: 


Vi (wt) = toc [cos (01) + cos (02) +--+ + cos (8p) ]sin (wt) (4) 


Where 01, 02,..., Ok are switching angles and must meet the following condition: 
0 <81 <02 <+ < Ok < 1/2 (5) 


In the SHE technique, the control of the fundamental voltage is carried out with the modulation index (M) 
which is the ratio of the peak value of the desired basic voltage (V ıp) to the total DC input voltage given in (6): 


aip. 
S eyes (6) 
The 3 harmonic components and multiples of 3 harmonics at phase-to-phase voltages can be 
disregarded in a balanced 3-phase system. As a result, the line voltage waveform's 5“, 7", and 11"... harmonic 
orders can be removed at the low switching frequency. In a 3-phase 11-level MLI, the following equations must 
be solved to eliminate undesired harmonics and generate fundamental harmonics of specified amplitude: 


V, = cos (81) + cos (02) + +--+ cos (65) = Mkr/4 

V; = cos(56,) + cos(50,) +++ cos(50;) = 0 

V, = cos(76,) + cos(762) + +--+ cos(76s) = 0 (7) 
Vi, = cos (116,) + cos (1103) + ---+ cos (110,) = 0 

Vız = cos (1304) + cos (136,) + = + cos (136,) = 0 


The fundamental waveform's amplitude is managed by the first SHE equation, as shown in (7). To 
remove the selected harmonics, other equations are set to zero. After that, these equations are correctly 
solved. The goal is to remove or significantly reduce the low-order harmonics (V5, V7, V11, and V13) while 
ensuring the fundamental frequency maintains a steady voltage value at the desired value. To find the optimal 
SHE equation solutions, an RDA optimization strategy is suggested in this study. 


4. IMPLEMENTATION OF RDA OPTIMIZATION TO SHE TECHNIQUE 

The RDA [25] was first reported; this was followed by the introduction of a new RDA optimization 
strategy [26]. The main source of inspiration for this meta-heuristic was the unusual mating pattern of Scottish 
red deers (RDs) during the breeding season. Like other meta-intuitive methods, RDA is initiated with a set of 
randomly selected RDs, where the best RDs in the population are considered the "male RDs," while the others 
are called "hinds." The male RD should roar first; they are split into two categories-single male deer (stags) and 
commanders based on the strength of the roar phase. The commanders engage the stags in a battle to claim their 
harems. Here, only the commanders can create the harems. The number of hinds in the harems of the leaders is 
directly related to their roar and their prowess in the battle phase. The commanders consequently mate a number 
of hinds in their harems but each commander is allowed a limited number of mating with hinds from various 
harems. The other male deers (stags) can also mate with the closest hinds irrespective of the size of the harem. 
Both exploration and concentration phases are also focused on in this step. The mating procedure by which RDs 
offspring are created is a significant stage in RDA optimization. This stage also offers new solutions to the 
optimization problem. The next algorithmic stage is finally created by considering the solutions that are 
evolutionarily weak [26]. The following steps are performed to solve the SHE problem with RDA: 
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a. Determination of the fitness function: the cost function for the SHE problem is defined as (8): 


f = ming; {In — (“2*) ; (ZI, cos (9;))| + (=. (x, cos cna))) | =0 (8) 


TT nT 


Where Vıp (H1) is the desired fundamental harmonic, hn; are the harmonics that need to be eliminated. The first 
term controls the fundamental voltage and the second term represents low-order harmonics in the cost function. 
b. Creation of the initial population: create the initial population according to (9): 


Value= f (Red Deer) = f (X1, X2, X3, ...,Xw,,,) (9) 


Where Xnvar Specifies the size of the array, while the values X1, X2, and X3. are the array's components. Five 

variables are defined (01, 62, 03, 84, and 05) since only one switching is made for each bridge in the eleven-level 

CHB-MLLI. The number of variables will be the number of positions assigned to each deer. To start the algorithm, 

a random initial population of Npop size is created, taking into account the fitness function. 

c. Choosing the best solutions: randomly generated switching angles are evaluated according to the fitness 
function. Some of the best solutions (depending on the value selected in Table 1) are identified and stored 
as male RD. 


Table 1. RDA optimization parameters 


Parameters Value Parameters Value 
Number of variables 5 Alfa (a) 1.0 
Iteration number 100 Beta (B) 0.6 
Number of the population 20 Gamma (y) 0.5 
Number of male RD 5 Roar 0.2 
Number of hinds 15 Fight 0.5 

Mating 0.8 


d. Selection of commanders and stags: select the commander and stags according to (10): 
Noom = round {y ` N male? (10) 
Nstag T N male ~ Ncom (11) 


Where Nmate is the number of all male deers, Ncom represents the number of commanders, and finally, Nstag is 
the number of stags (single male deers). 
e. Fight: each commander randomly fights the stags. The fighting process is performed using (12): 


_ (Com+Stag) 


New, = ———— + b, x (((UB — LB)* b2) + LB) 


ce (12) 
=o -h x (((UB — LB)* bz) + LB) 


New, 


Where UB is the upper bound of the search space while LB is the lower bound. The fighting step is a random 

process and as such, the selection of bı and b2 is done randomly between 0 and 1 in a uniformly distributed 

manner. A commander can approach a stag or vice versa in the resolution area. Any of the four solutions 

[single male RD (Stag), commander (Com), New1 (new status1), and New2 (new status2)] with better FF 

becomes the new commander. 

f. Create harems: a harem is made up of a male commander and a group of captured female deer (hinds). 
The fighting capability of the commander determines the number of female deer in the harems. 

g. Mating: mating can be done in three ways: i) mating each commander with females from his harem; ii) 
mating any commander with hinds from another harem; and iii) mating the stag with the nearest hind, 
irrespective of the harem restrictions. 

h. Creation of offspring: during the mating phase, offspring are produced in all three cases according to (13). 
The Stag replaces Com during the third case in the mating phase: 


(Com+Hind) 
2 


ofss = + (UB — LB) xc (13) 


Where c is randomly generated in a uniformly distributed manner between 0 and 1. 


i. Selection of the new generation: the next generation is chosen using two alternative methods. The first 
method involves keeping the best male RDs (Stag and Com) (a given proportion of the optimal solutions). 
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The second technique requires the selection of the hinds for the next generation utilizing a fitness 
competition or roulette wheel mechanism; the selection is done from among all the hinds and produced 
offspring during the mating stage. The selected male and female RDs produce the next generation. 

j. Stop criterion: the number of iterations, the best-found solution, or a given time interval can all be 
considered as the stopping condition. The number of iterations is used as the stopping condition for 
RDA-based SHE-PWM. Repeat steps c through j until the stopping condition is satisfied. 

Details of how RDA optimization can be used with SHE equations have been provided by 

[23], [24]; readers are referred to the relevant articles for more explicit information. RDA optimization 

parameters are adjusted to solve SHE equations in 11-level three-phase CHB-MLI. Figure 2(a) shows the 

optimum switching angles that the RDA optimization found in the range of modulation index of 0.1 to 1.0 to 
solve the SHE equation. Figure 2(b) shows the optimum switching angles that the proposed method finds to 
provide the OMTHD value. Figures 3(a) and (b) show total harmonic distortion (THD) and THD, values within 
the modulation index range of 0.1 to 1.0 for SHE and OMTHD, respectively. Evidently, the proposed RDA- 
based method finds suitable solutions for SHE and OMTHD in the modulation index range of 0.4 to 1.0. 


RDA Optimization (SHE) 
T T T 


RDA Optimization (OMTHD) 
r T T T r 


Switching angles degrees (8) 
Switching angles degrees (0) 


1 f f f 1 f 1 J 1 1 f 1 f f 
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 07 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Modulation Index (M) (a) Modulation Index (M) (b) 
(a) (b) 


Figure 2. Switching angles according to modulation index (a) for SHE and (b) for OMTHD 
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Figure 3. THD% and THD.% according to modulation index (a) for SHE and (b) for OMTHD 


Two types of THD are computed in this work; the first is the overall THD; the second is the THD. 
value that stands for the total of harmonics until the eliminated maximum harmonic orders. In the 11-level 
three-phase CHB-MLI, THD, and THD, can be calculated as given in (14) and (15), respectively. The upper 
limit for THD is infinite, but THD, is calculated up to the minimized harmonic orders. 


\v2+v24+V2, + 
THD = +——_—_ 


Iva] 


(14) 


2 2 2 2 
v +V AVE H+ Ve a 


THD, = , (for odd k.) (15) 


Iva] 


5. RESULT 
To prove the effectiveness of the proposed algorithm, the calculated switching angles were applied 
to the three-phase 11-level CHB-MLI. Each H-bridge module input voltage is selected as 50 V (Vac=50 V). 
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These results show that the selected harmonics are minimized almost to zero and the fundamental harmonic 
is kept very close to the desired value. The switching angles in the unit modulation index to be minimized of 
selected harmonics are estimated as 01=6.75°, 02=18.98°, @3=27.53°, 04,=45.49°, and 05=62.40°. The 
switching angles in the unit modulation index for the best OMTHD value are obtained as 01=8.02°, 02=19.2°, 
83=27.82°, 04=46.45°, and 65=61.16°. For the best SHE solution, the inverter output phase voltage waveform 
is given in Figure 4(a). The frequency spectra of line voltage for THD and THD, are depicted in 
Figures 4(b) and 4(c), respectively. Similarly, the output phase voltage and the frequency spectra of line 
voltage for THD and THD, are given in Figures 5(a)-(c), respectively. 

For SHE and OMTHD, the fundamental phase voltage is normally V;=250 V for M=1.0 from (6). 
The switching angles, THD, and THD. values calculated for the NR, LSHADE, WOA, PSO, and the 
proposed RDA-based technique are given in Table 2. As can be seen, the proposed RDA-based SHE-PWM 
technique (RDA-SHE) eliminates selected harmonics better than another optimization techniques. The 
smallest THD, value proves that the selected harmonics are efficiently minimized. 

The table suggests that when the calculated switching angles for the RDA-SHE are applied, the 
THD,.% value is obtained as 0.04% and the THD% value as 4.50%. By using the suggested RDA method, the 
amount of intensification can be adjusted between the harmonics desired to be eliminated and the THD value. 
If prompted to obtain the best THD value, the intensification of the THD can be increased in the fitness 
function. In this case, the THDe% and THD% values are obtained as 1.29 % and 3.96 %, respectively. From 
the perspective of minimum THD, the proposed RDA has given the best performance, but the harmonic value 
is higher than the calculated value for the SHE because optimization (RDA-OMTHD) is intensified to find a 
solution for the minimum THD value. 
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Figure 4. For SHE solution (a) phase voltage waveform, (b) THD value of line voltage, and (c) THD, value 
of line voltage 
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Figure 5. For OMTHD solution (a) phase voltage waveform, (b) THD value of line voltage, and (c) THD. 
value of line voltage 


Table 2. Comparison of THD% and THD.% values of different algorithms 


Switching angles (degree) 


Optimization techniques Modulation index (M) el @2 o3 o4 o5 THD% THD.% 
NR method [16] 1.00 6.56 18.93 27.17 45.12 61.87 4.35 0.26 
PSO [18] 1.00 6.72 19.00 27.39 45.10 62.16 4.45 0.17 
LSHADE [19] 1.00 5.49 16.68 28.58 42.06 59.46 4.67 1.83 
WOA [20] 1.00 6.60 18.80 27.30 45.10 62.20 4.53 0.09 
Proposed RDA-SHE 1.00 6.75 18.98 27.53 45.49 62.40 4.50 0.04 
Proposed RDA-OMTHD 1.00 8.02 19.28 27.82 46.45 61.16 3.96 1.29 
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6. CONCLUSION 

In this paper, an RDA optimization is proposed for the elimination of selected harmonic orders in 
three-phase 11-level CHB-MLI fed by symmetric sources. Suitable switching angles for the modulation 
index values between 0.1 and 1 for the suggested RDA optimization were quickly calculated. The results 
show that the RDA method can successfully attenuate low-order harmonics while maintaining the 
fundamental output voltage. Furthermore, the proposed method can be used effectively to achieve the optimal 
THD value. Consequently, the proposed optimization method produces better solutions within the specified 
modulation index range compared to other existing optimization methods. 
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